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Renormalized viscosity, renormalized resistivity, and various energy fluxes are calculated for he- 
lical magnetohydrodynamics using perturbative field theory. The calculation is to first-order in 
perturbation. Kinetic and magnetic helicities do not affect the renormalized parameters, but they 
induce an inverse cascade of magnetic energy. The sources for the the large-scale magnetic field 
have been shown to be (1) energy fiux from large-scale velocity field to large-scale magnetic field 
arising due to nonhelical interactions, and (2) inverse energy flux of magnetic energy caused by 
^^ ' helical interactions. Based on our flux results, a premitive model for galactic dynamo has been 

f^ ' constructed. Our calculations yields dynamo time-scale for a typical galaxy to be of the order of 

("^ , 10* years. Our field-theoretic calculations also reveal that the fiux of magnetic helicity is backward, 

rvQ ' consistent with the earlier observations based on absolute equilibrium theory. 

Q^ PACS numbers: PACS numbers: 47.27.Gs, 52.35.Ra, 91.25.Cw 
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lO ' I- INTRODUCTION 

(N. 

,__, , Generation of magnetic field in plasma, usually referred to as "dynamo" , is one of the prominent and unsolved 

^**^ ■ problems in physics and astrophysics. It is known that the magnetic field of galaxies, the Sun, and the Earth are 

r -v , neither due to some permanent magnet nor due to any remnants of the past, but it is generated by the nonlinear 

^^ ' processes of plasma motion ('1, 2] and references therein). However, a solid quantitative understanding is lacking in 

fi , this area inspite of various attempts for more than half century. There are various aspects in this problem, and we 

1— H . address energy transfer issues in this paper and in paper 1(3]) using field-theoretic methods in a somewhat idealized 

, " , ' environment, homogeneous and isotropic flows. In paper I we show that the nonhelical part of magnetohydrodynamic 

(MHD) interaction causes energy transfer from large-scale (LS) velocity field to large-scale (LS) magnetic field. In a 

^^ ' typical dynamo environment however, the helical interactions cause an additional energy cascade of magnetic energy 

K*" , from small scales (SS) to large scales; the field-theoretic calculation of helical contribution to the energy flux is 

^^ presented in this paper. Both helical and nonhelical factors contribute to the magnetic energy growth. 

In the problem of magnetic field generation, it is required that the LS magnetic field is maintained at all time. 
There are several exact results in this area, e.g., dynamo does not exist in two dimensions as well as in axisymmetric 
flows [l|. In the past, several dynamos, e.g., rotor dynamo, 2-sphere dynamo etc., have been constructed [l|, however, 
mean- field electrodynamics developed by Steenbeck et al. p] (also see Krause and Radler ^) paved a way for 
f^ ' practical calculations in astrophysical and terrestrial dynamo. This formalism also provided insights into the physical 
mechanism of dynamo, mainly that kinetic helicity Hk = l/2((u • Jl), where u and Jl are the velocity and vorticity 
fields respectively) plays an important role in the amplification of the magnetic field. The amplification parameter 
(— I , Uu was found to be {t/3)Hk, where r is the velocity de-correlation time. 

Mean-field electrodynamics of Steenbeck et al. U is a kinematic theory. Here it is assumed that the velocity field is 
a known function, which is unaffected by the generated magnetic field. The later models which take into account the 
back reaction of the magnetic field to the velocity field are called dynamic models. One of the first dynamic model 
$H ' is due to Pouquet et al. J5| where they incorporated the feedback, and proposed that the modified a is proportional 
5^ , to residual helicity, Hk — Hj, where Hj is the current helicity define as 1/2 (b • Vb). Gruzinov and Diamond Jg 
proposed a quenching mechanism where the amplification parameter was modified to 

(1) 
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where B is the LS magnetic energy, U is the LS velocity field, and Rm is the magnetic Reynolds number. Recently 
Field, Blackman, and Chou 7], and Chou [g obtained a general expression for dynamic a coefficient as a function 
of Reynolds number and magnetic Prandtl number. Basu and Bhatthacharya [3 and Basu 10] have attempted to 
compute the dynamo coefficients a and (3 using field-theoretic techniques. 

In a recent development Brandenburg ,20] investigated dynamo problem in isotropic and helical MHD. The the 
system is forced with kinetic helicity, they find magnetic energy transfer to large-scales. They indentify this mecha- 
nism, named as nonlinear alpha-effect, for magnetic energy growth. Magnetic helicity play an important role in this 
mechanism; we will discuss these issues in later part of the paper. 

There are many numerical simulations of dynamo in various geometries. In simulations MHD equations are numer- 
ically solved with appropriate boundary conditions. In this paper we will only refer to the simulations performed for 



a periodic box; this is to avoid the complications of spherical geometry (e.g., effects of Coriolis force etc.). Pouquet et 
al. 13 numerically integrated the MHD equations on the basis of Eddy-damped-quasi-normal-Markovian (EDQNM) 
approximation. When kinetic energy (KE=l/2 (u • u)) and Hk were injected near a wavenumber band, both mag- 
netic energy (ME=l/2 (b • b)) and absolute value of magnetic hclicity \Hm\ were found to increase. Magnetic helicity 
Hm is defined as 1/2 (a • b), where a is vector potential. Note however that one part of the time-scale used in Pou- 
quet et al.'s calculation is based on Alfven relaxation time. This assumption is suspect in view of current theoretical 
[llL [l3 . [13I [Q. flSJ I and numerical results p^[l7llla |. which favour k~^/^ (Kolmogorov's) energy spectrum over Kraich- 
nan's k~^''^ spectrum for MHD turbulence. In our present paper, the nonlinear time-scale is based on Kolmogorov's 
energy spectrum. The direct numerical simulation (DNS) of Pouquet and Patterson [13 yielded a similar result. The 
results of Brandenburg [23| discussed above are also obtained numerically. Frick and Sokoloff [21| studied the shell 
model of turbulence, and showed that magnetic helicity suppresses turbulent cascade. 

In another development Dar et al. |22!| numerically calculated various energy fluxes in two-dimensional MHD and 
showed that there is significant energy transfer from LS velocity field to LS magnetic field; they claimed that the above 
flux is one of the main contributor to the amplification of large-scale ME. Note that the ratio ME/KE grows in both 
helical and nonhelical MHD as well as in many decaying as well as forced simulations (see for example, |22, 23] and 
references therein). Hence, helicity is not a necessary requirement for the generation of magnetic energy. Regarding 
flux of magnetic helicity, Pouquet et al. [3 and Pouquet and Patterson [13 argue that it is in the inverse direction 
(from SS to LS). We find a similar magnetic helicity flux in our theoretical calculation. 

In one of the recent theoretical development Kulsrud and Anderson 24] derived and solved the kinetic equation 
for the growth of galactic magnetic field. They argued that the dynamo time-scale is much large than the growth 
time-scale of turbulent modes. Hence, the buildup of SS turbulent ME dominates the slow growth of LS ME, thus 
making the mean dynamo theory invalid. 

In paper I and the present paper, the energy fluxes of MHD turbulence are computed using perturbative field theory. 
The calculation is to first order in perturbation. Here the viscosity and resistivity were taken from renormalization 
calculation of Verma [l4lll5l |. and the energy spectrum were taken to Kolmogorov-like (k~^/'^). Kolmogorov's spectrum 
for MHD turbulence is supported by recent theoretical [l^, [l3; HM Hji 1^ ^^nd numerical results jla, [l3i US • In paper 
I we show that in a kinetically forced nonhelical MHD, the energy transfer from LS velocity field to LS magnetic field 
is one of the dominant transfers. It is also shown in paper I that the above energy flux into LS magnetic modes is 
independent of the nature of LS forcing. In this paper we generalize the field-theoretic calculation of Verma Q, 13 
to helical MHD. As discussed in paper I, we assume that the turbulence is homogeneous and isotropic, and that the 
mean magnetic field is absent. The absence of mean magnetic field is a reasonable assumption for the initial stages 
of dynamo evolution. This assumption is to ensure that the turbulence is isotropic. In addition we take cross helicity 
(u • b) to be zero to simplify the calculation. We examine the fluxes in the presence of magnetic and kinetic helicities. 
We also investigate the flux of magnetic helicity. 

We have constructed a simple dynamo model using the theoretically-calculated energy flux into the large-scale 
magnetic field. In this model the ME grows exponentially in the initial stage, and the growth time-scale is of the 
order of eddy turnover time. It shows that dynamo action is possible in galactic dynamo. 

The outline of the paper is as follows: in section 2, we carry out the perturbative calculation of renormalized 
viscosity and resistivity, as well as fluxes of energy and magnetic helicity. In section 3, we construct a dynamic 
galactic dynamo based on our energy flux results. In this section we also compare our results with the findings of 
earlier researchers. Section 4 contains conclusions. 

II. FIELD-THEORETIC CALCULATION OF HELICAL MAGNETOHYDRODYNAMICS 

The incompressible MHD equation in Fourier space is given by 

[-iuj + vk^) u,(k) = --P.;^^{\i.) dp[uj{p)u„^{k-p)-bj{p)b,n{k~p)] (2) 

{-iLu + r^e)bdk) = -iPr^^ik) J dp[uj{p)bm{k-p)] (3) 

hu,{k) = (4) 

kAik) = (5) 

where u and b are the velocity and magnetic field fluctuations respectively, v and 77 are the viscosity and the resistivity 
respectively, and 

i^+„(k) = k,PU^) + k^P^.ik); (6) 
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Note that we are working in three space dimension and also with zero mean magnetic field. 

Some of the definitions regarding kinetic and magnetic helicities are in order. Throughout this paper, a denotes the 
vector potential (b = V x a) , and fJ denotes the vorticity (fi = V x u) . The spectrum of helicity, Hm (k) , is defined 
using the equal-time correlation function (ai(k, t)6j(k, t)) (the angular brackets denote ensemble average), 

(a,(k,i)6,(k',i)) =>y(k)i/M(k)(27r)35(k + k') (10) 

The factor -Pij(k) appears due to the constraints V • a = V • b = 0. Using this correlation function we derive the 
following relationship: 

Hm = ^(a(x).b(x)) 

' (a(k).b(k')) 



(27r 
The one-dimensional magnetic helicity HM{k) is defined using 
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-^HM{k) (11) 



Therefore, 



Using V X a = b we can easily derive that 



dk 

HM{k)dk = / j^HMik) (12) 



HM{k) = ^HM{k) (13) 



a(k) = ^k X b(k) (14) 



which leads to 

(6,(k)fe,(k')> = [Py (k)C''^(k) - ie,,ikiHM{k)] {2TTfd{k + k') (15) 

where C^^ is the b — b correlation function. 

A similar analysis for kinetic helicity shows that 

{v,{k)n,{k')) 1^' P,,{k)HK{k)Sik + k') (16) 

2^ 'J (27r) 



i/K = - (u • O) = / j,^-^HKik) (17) 



and 

HK{k) 



{u,{k)u,{k')) 



P„(k)C""(k)-ze,,,fcr 



{2nf6{k + k') (18) 



where C"" is the u~ u correlation function. 

An important point to note is that magnetic helicity is conserved in MHD when v — rj = Q. One of the consequences 
of this conservation law is the emergence of k~^ spectrum at smaller wavenumbers 25]. The kinetic helicity is conserved 
in fluid turbulence, but not in MHD turbulence. 

In the following subsection we will calculate the renormalized viscosity and resistivity for helical MHD. 



A. Calculation of renormalized parameters 

Recently Verma [l^ [l^ has calculated the renormalized viscosity and resistivity for MHD turbulence in absence of 
kinetic and magnetic helicity. It will be shown below that the presence of both kinetic and magnetic helicities does 
not alter the renormalized viscosity and resistivity calculated for nonhelical MHD. 

In the RG procedure the wavenumber range {kpf, kg) is divided logarithmically into N shells. Then the elimination 
of the first shell k^ — (fci,/co) is carried out and the modified MHD equation for A:< = {kpf,ki) is obtained. This 
process is continued for higher shells. The shell elimination is performed by ensemble averaging over fc^ modes 
[Tj, [13, 123 . It is assumed that uf{k) and bf{k) have gaussian distributions with zero mean, while uf(k) and bfik) 
are unaffected by the averaging process. In addition it is also assumed that 
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(19) 
(20) 



Note that u — b correlation has been taken to be zero in our calculation. 

We apply first-order perturbation theory to compute the renormalized parameters. After elimination of n shells, we 
obtain the following equations for the renormalized viscosity vtn) ^-nd renormalized resistivity rirn\ (for details refer 
to Verma ^). 



-iuj + V(n)k'^ + (5i^(„)fc^) uf{k) 



^iuj + ?7(„)fc^ + 5-q{.a)k^) bf{k) 



'^Pijv^O^) j dp[uf{p)u<{k-p) 
-b<{p)b<{k-p)] 
-*/^_(k) fdp[u<{p)b<ik-p)] 



where 
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S{k,p,q)C^^{q) + S'{k,p,q)HK{q) 
V(n){p)P^ + V(n){q)q'^ 

S,{k,p,q)C^\q) + S'^{k,p,q)HM{q) 



dp 

'p+q=k (27r)3 



'n(n){p)p'^ + %n){q)q^ 

Ss{k,p,q)C'\q) + S's{k,p,q)HM{q) 

v(n){p)p^ + V(n){q)q'^ 

, Ss,{k,p,q)C^^{q) + S'^(k,p,q)Hk{q) 



V(n){p)p'^ + l'(n) (9)9^ 



The quantities 5,; and S[ are as follows: 

S{k,p,q) = P,+,Jk)P+^,ip)P,a{q)^2kp{z' 



xy 



Seik,p, q) = P^j^{k)p-,,Sp)PMi) = -2fcp^ (1 - y') 

S^{k,p,q) = Pr^^{k)P+^,{p)P^MP^b{k) ^ S^{p,k,q) 

S,{k,p,q) - Pr^.,^{k)p-^^,{p)P,a{q)P^b{k)^'^kp{z + xy) 

S'{k,p,q) = P,+Jk)P+Jp)e,aiqi=0 

S'e{k,p,q) = PajrrXk)P,nba(P)<^Jaiqi ^ 

S',ik,p,q) = PiJ^Jk)P+,ip)e^alqlP^bik)=0 

Sgik,p,q) = Pij,nik)PmabiP)^J'^l<llP^bik) ^ 



(21) 
(22) 



(23) 



(24) 



(25) 
(26) 
(27) 
(28) 
(29) 
(30) 

(31) 

(32) 



Since Siy and Sr] are proper scalars and Hm.k are pseudo scalars, S^{k,p,q) will be pseudo scalars. In addition, 
S[{k,p, q) are also linear in k,p and q. This implies that S[{k^p, q) must be proportional to q • (k x p), which will be 
zero because k = p + q. Hence all S[{k,p, q) turn out to be zero, as a consequence the presence of helicities does not 
alter the already calculated (5(;/, ??)(„) (fc) by Verma [ij,[l3- Zhou [23 arrived at a similar conclusion while calculating 
the renormalized viscosity for helical fluid turbulence. 



Verma [ij, [T^ obtained a self-consistent solution of the renormalized parameters using Kolmogorov's spectrum. 
Since the helicities do not alter the renormalized parameters, we arrive at the same formula for renormalized viscosity 
and resisitivity as Verma |lj,[i3, that is, 

( w,, _ / (if")'/'ni/3fc-4/3(^*,^*) forfc>fc„ 

^^^^m - I (if«)l/2nl/3fc,-4/3(^*^^*) for k < K ^ ' 

where 11 is the total energy flux, X" is the Kolmogorov's const ant, and v* and rj* are the renormalized parameters. 
The value of these renormalized parameters have been listed in |l4l Il5l | . 

The present calculation has been carried out up to first order. The probability distribution of velocity is gaussian, 
while that of velocity difference is nongaussian. The nongaussian behaviour of velocity difference has significant 
effects specially on higher order structure functions, which are not properly accounted for by first order calculations. 
Yet, the first order calculation of renormalized viscosity yields results very close to those obtained in experiments 
and numerical simulations (see e.g., |2a. l2q'). For the above reason, we have stuck to the first-order field-theoretic 
calculation in the present paper. 

In the next subsection we will calculate the energy and helicity fluxes using the field theoretic technique. 

B. Calculation of energy and helicity fluxes 

In paper I we have analytically calculated energy fluxes in the absence of magnetic and kinetic helicities. In this 
subsection we will generalize that calculation for helical MHD. Refer to paper I for the energy evolution equations 
and other basic formulas. 

As discussed in paper I, the energy flux from inside of the X-sphere (X <) to outside of the F-sphere (F >) is 



^y>(ko) = ,,^,3.n. ^ / tSs / 7S3(^^^(k'|p|q)) (34) 



1 r dk' r _rfp_ 

(27r)3<5(k' + p + q) A,>,„ (2^)3 J^^,^^ (2^)3 

where X and Y stand for u or b, and 5'(k'|p|q) is energy transfer from mode p of X field to mode k of F field, with 
mode q acting as a mediator. The detailed expressions for (S'^'^(k'|p|q)) are given in Paper I. 

We calculate the above fluxes analytically to the leading order in perturbation series using the same procedure as 
in paper I. Some additional terms appear in {S{k'\p\q)) due to the presence of helicity. The detailed expressions are 
given in Appendix A. 

A formula for the magnetic helicity flux can be derived in a similar manner. From Eqs. I|2I3|I we can easily obtain 
the equation for the evolution of magnetic helicity, which is 



dt 2(27r)3,5(k -h k') 



^,,, , aa(k') ,., , ab(k') 
b*(k)-^i-^+a*(k) 



dt ^ ' at 



(35) 



where 



- 2(2.)3/(k + kO [^""(>^'IPl^) + ^"-(k'lqlp)] (36) 

5^-(k'|p|q) = i3?[b(k')-(v(p)xb(q))] 

+ ^9 [k' . b(q)a(k') . u(p) - k' . u(q)a(k') • b(p)] (37) 

Here 5R() and 3() stand for real and imaginary part of the arguments, respectively. The quantity S'^"^(k'|q|p) can 
be obtained from the above expression by interchanging p and q. After some algebraic manipulation it can be shown 
that 

5^'^(k'|p|q) + 5^^-^(k'|q|p) + 5^^'(p|k'|q) 
+5^"(p|q|k') + 5^*^(q|k'|p) + 5^^'(q|p|k') = (38) 

It shows that the "detailed conservation of magnetic helicity" holds in a triad interaction (when v = r] = 0) |23 • 
From the above, the transfer rate of magnetic helicity from a wavenumber sphere of radius fco is 



We again compute ^^^^^(k'lplq)) to first order in perturbation. The detailed expressions are given in Appendix B. 
The expressions in the Appendices involve Green's functions and correlation functions. The expressions for these 
functions are taken from self-consistent calculations (see e.g., Verma [l^)- For G(fc, t — i') of the formulas (|A1IA4IB3|I . 
we substitute 

Q(uu,bb) (^^ ^ _ i') = 5,(i _ ^') exp [- (i^(fc), ?7(fc)) fc2(i - t')] (40) 

where {v{k),ri{k)) are given by Eq. (|33|l . and 9{t — t') is the step function. We assume the relaxation time-scale 
for C'^'^{k,t,t') and HK{k,t,t') to be {L'{k)k'^)-\ while that of C^^{k,t,t') and HM{k,t,t') to be {ri{k)k'^)-'^ . The 
spectrum C^"'^'^''\k,t,t) are written in terms of one-dimensional energy spectra E'^^'^'> as 

^iuu,bb) (^^ ^^ ^) ^ I^£;(«,fc) (41) 

In presence of magnetic helicity, the calculations based on absolute equilibrium theories suggest that the energy 
cascades forward, and the magnetic helicity cascades backward Q^ . In this paper we have not considered the inverse 
cascade reg ion of magnetic helicity. We take Kolmogorov's spectrum for energy based on recent numerical simulations 
[la.ll7l ll q | and theoretical calculations pTl Il3. Il3lll5l | (ignoring the intermittency corrections). Hence, the spectrum 
of £;("''') can be taken as 

£;"(/c) = X"n2/3fc-5/3 (42) 

E\k) = E^JTA (43) 

where 11 is total energy flux. 

The helicities are written in terms of energy spectra as 

llK{k) - TKkE^ik) (44) 

HM{k) = m— P (45) 

We are calculating energy fluxes for the inertial-range wavenumbcrs where the same powerlaw is valid for all energy 
spectrum. Therefore, the ratios r^,rM, and rx can be treated as constants. 

We substitute the above forms for the correlation and Green's functions [Eqs. (|40I43|) ] in the expressions for 
(S^^ {k'\p\q)) and (S^" (k'lplq)) given in the Appendices. These 5"s are substituted in the flux formulas (Eqs. |^ 
139]'). We make the following change of variable: 

, ^0 ^0 ko , , 

fc = — ;p= — v:q~ — w (4d) 

u u u 

These operations yield the following nondimensional form of the equation in the —5/3 region (for details, refer to j^j). 

n?^(fco 



n(fco 



= (if")3/2 



ngA/(fco) ^ J_CJ^u\3/2 



1 pl+v 

dvh\{l/v) I dw{vw) sin aFyy 

Jl-ii 

1 ri pi+v 



2./0 



dv{l ^ v) / dw{vw) sin aFn,, 



(47) 
(48) 



where the integrands {Fyy ^Fhm) are function of u, w, u* , rj* , rA,rK and tm I^- 

We compute the term in the square brackets, /y> , using the similar procedure as that of Verma 3]. The flux ratios 
Ily ^ /n can be written in terms of integrals ly^ , which have been computed numerically. Table Q] contains their 
values for ta ~ I and rA = 5000. The constant K'^ is calculated using the fact that the total energy flux 11 is sum 
of all n^<. For parameters (r^ = 5000, r^ = 0.1, a/ = -0.1), i^" = 1.53, while for {ta = I.tk ^ 0.1, tm = -0.1), 
-fC" — 0.78. After this the energy flux ratios Ily^/II can be calculated. These ratios for some of the specific values 
of rA, rx and tm are listed in Table Ull The flrst and second terms of Ily/II entries are nonhelical and helical 
components respectively. 

An observation of the results shows some interesting patterns. The energy flux can be split into two parts: helical 
(dependent on vk and/or tm) and nonhelical (independent of helicity). The nonhelical part of all the fluxes except 
nt> (n^>nonfteiicai < ^ ^^^ ^A > 1) is always positive. As a consequence, in nonhehcal channel, ME cascades from LS 
to SS. Also, since 11^^ > 0, LS kinetic energy feeds the LS magnetic energy. The fluxes of nonhelical MHD has been 
discussed in great detail in paper I. 



The sign of ^u>heiicai is always negative, i.e., kinetic helicity reduces the KE flux. But the sign of heUcal component 
of other energy fluxes depends quite cruciaUy on the sign of helicities. From the entries of Table P we see that 

^\b>,u»heHcal = -^^I + ^^MTK (49) 

where a and h are positive constants. If tmtk < 0, the energy flux to LS magnetic fleld due to both the terms in the 
right-hand-side of the above equation is positive. Earher EDQNM |^ and numerical simulations [23 with forcing of 
KE and Hj^ typicaUy have t^tm < 0. Hence, we can claim that helicity typicaUy induces an inverse energy cascade 
via n|j^ and n(^>. These fluxes wiU enhance the large-scale magnetic field. 

From the entries of Table UTI we can infer that the for small and moderate rK and rjv/, the inverse energy cascade into 

large-scale magnetic field is less than the forward nonhelical energy flux H^-,. While for helical MHD (r^, tm ^1), 
the inverse helical cascade dominates the nonhelical magnetic-to-magnetic energy cascade. 

The flux ratio Hhm /^ can be written in terms the integrals of Eqs. H48I48|I using the same procedure as done for 
energy flux ratios. The numerical values of the integrals are shown in Tables 1 and 2. Clearly, 

^Hm = -drM + erK (50) 

where d and e are positive constants. Note however that contribution of Hm dominates that of Hk- Clearly the sign 
of ^Hm is the same as that of Hk but negative of Hm- From the above equation we observe that positive Hm yields 
a negative contribution to ^Hm- Hence, for positive Hm, the magnetic helicity cascade is backward. This result is 
in agreement with Frisch et al.'s |23 argument in which they predict an inverse cascade of magnetic helicity. Our 
theoretical result on inverse cascade of Hm is also in agreement with the results derived using EDQNM calculation 
[3 and numerical simulations [l9| . 

When we force the system with positive kinetic helicity [tk > 0) , Eq. (|50|l indicates a forward cascade of magnetic 
energy. This effect could be the reason for the observe production of positive magnetic helicity at small scale by 
Brandenburg 20i| Because of magnetic helicity conservation, he also finds generation of negative magnetic helicity at 
large-scales. Now, positive kinetic helicity and negative magnetic helicity at large-scales may yield an inverse cascade 
of magnetic energy (see Eq. |49|l . This could be the crude reason for the growth of magnetic energy in the simulations 
of Brandenburg [23 • 

In paper I we calculated H^^ for nonhelical MHD using steady-state condition. 

TTM< _ TTfc< -I- TT''< (^^'] 

^^fc< '^'^bynonhelical ^ '^''uynonhelical- \^^J 

The above calculation for helical MHD is not straight forward because magnetic energy at large-scale could increase 
with time, and steady state may not be achievable for all possible parameters of helical MHD. Brandenburg ^^ 
observes the dynamic evolution of large-scale magnetic energy in his simulations. To simplify the calculation, we 
assume steady-state condition for helical MHD as well, and calculate various parameters. For some set of highly 
helical MHD, we get negative energy fiux. 

In the following section, we will construct a dynamic dynamo model for galaxies using our flux results. 

III. DYNAMO VIA ENERGY AND MAGNETIC HELICITY FLUXES 

In the above calculation we have assumed that the turbulence is homogeneous, isotropic, and steady. The assumption 
of homogeneity and isotropy can be assumed to hold in galaxies in the early phases of evolution before large structures 
appear. The assumption that the mean magnetic fleld of galaxy is rather small is valid in the beginning of the galactic 
evolution. Therefore, we apply the flux obtained from our calculations to estimate the growth of magnetic energy in 
galaxies. 

During the early phase of galactic evolution, only the large-scales (LS) contain the kinetic and magnetic energies. 
The fields at these scales interact with each other, but the small-scale spectrum is far from steady (not enough time). 
The interactions of LS velocity field and the LS seed magnetic field increase the LS seed magnetic field E^{t) till the 
steady-state is reached. In absence of helicity, the source of energy for the large-scale magnetic field is H^^ [Eq. (|51|) ]. 
When helicity is present, there are several other sources as discussed in section 2 of this paper. Since the forcing of 
helicities is effective at LS, it is reasonable to assume that the helical part of H^^ and H|^^ will also aid to the increase 
in LS ME. Hence, 






nfc< + ^bKhelical + ^b<helical (52) 
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We assume a quasi-steady approximation for the early evolution of magnetic field. In many quasi-steady situations 
(slowly decaying or growing), steady-state results are usually applied. This approximation works very well for many 
practical problems. We make this assumption in this paper, and substitute the theoretically calculated energy fluxes 
calculated in Section II of the paper to the above equation. 

Since the ME starts with a small value (large ta limit), all the fluxes appearing in Eq. H52I) are proportional to r~^ 
[cf. Eqs. IA2IA4II ]. i.e., 

n^fx + '^b<helical + '^b<helical = ~£u '°'^) 

where i?" is the LS KE, and c is the constant of proportionality, which depends on the values of helicities. Both 
Kinetic and magnetic helicities are difficult to ascertain for a galaxy due to lack of observations. We take tm and r^ 
to be of the order of 0.1, with tm being negative. The choice of negative tm is motivated by the results of EDQNM 
calculation 5] and numerical simulations |23|. With this value of tm and r^, c ~ 0.84 for E{k) ex fc~^/^ regime, and 
c ss 1.3 for Eik) oc k~^ regime. Since both the values of constant c is approximately equal and close to 1.0, we take 
c — 1.0 for our calculation. Hence, 

1 dE^ E' (54) 



n dt £;« 

Using E'^ = K^lP/^L^/"^, where L is the large length-scale of the system, we obtain 

1 dE'' 1 



^JWE^ dt LiK^fl"^ 
We assume that E" does not change appreciably in the early phase. Therefore, 



(55) 



E\t)^E\0)e^^yj^^^^^t^ (56) 

Hence, the ME grows exponentially in the early periods, and the time-scale of growth is of the order oi LiK"^)^ '"^ j \fE^ ^ 
which is the eddy turnover time .3J. Taking L « lO^^fcm and \J E"^^ « 10km/ sec, we obtain the growth time-scale to be 
lO^^sec or 3 X 10^ years, which is in the expected range [2J|. Hence, we have constructed a nonlinear and dynamically 
consistent galactic dynamo based on the energy fluxes. In this model the ME grows exponentially, and the growth 
time-scale is reasonable .24]. 

The helical and nonhelical contribution to the fluxes for ta — 5000, rx — 0.1, tm — —0.1 is shown in Table ITll The 
flux ratios shown in the table do not change appreciably as long as ta > 100 or so. The three fluxes responsible for the 
growth of LS ME are U]^< /U ^ 2.6 x lO'^ (nonhelical), H^<^,,_,/n « -4.1 x lO'^, and H^< ,,,,,,,,/H « -4.0 x 10"^ 
The ratio of nonhelical to helical contribution is 2.6/0.81 « 3.2. Hence, the nonhelical contribution is signiflcant, if 
not more, than the contribution from the helical part for the LS ME amplification. Note that in the earlier papers 
on dynamo, the helical part is strongly emphasized. 

Kulsrad and Anderson (KA) |2J| performed an important mean field dynamo calculation of galactic dynamo for 
large Prandtl numbers. Some of the salient features are follows. In KA's kinematic dynamo calculation the growth 
rate of ME is 7 = (1/3) / k^U{k)dk where 

Cn'l2/3 1,-5/3 

KA estimate 7 « kmaxO^y^"^ ~ lO^^yr^^. From this result KA conclude that the kinematic theory predicts a 
extremely rapid growth of SS ME. The SS noisy magnetic field thus generated will dominate the mean magnetic 
field that grows at a considerably slower rate (dynamo growth time « 3 x lO^yr). Therefore, it is claimed that the 
kinematic assumption of the mean dynamo theory is invalid, and it is difficult to build up galactic magnetic field 
from a very weak seed field using dynamo action. KA's estimate of growth time-scale is equal to the eddy turnover 
time of smallest eddies {k^\^). Hence, as pointed out by KA, kinematic assumption is invalid for galactic dynamo, 
and one has to resort to a dynamical model. Brandenburg's numerical results 20] are not quite consistent with KA's 
results. For example, Brandenburg finds (1) growth of magnetic energy even for large magnetic Prandtl number; (2) 
the growth time-scale for magnetic energy is of the order of L'^/t], where L is the large length-scale, and 77 is resistivity. 
Our results are valid for Prandtl number close to 1. Therefore, they can not be compared with KA's calculations. 
It is interesting to see however that our crude estimate of time scale is one eddy turn-over time. To get a better 
picture, we need to construct a more solid model. 



In our model the magnetic energy growth is due to the fluxes n|^< + ^b<heiicai + ^b<heiicai- ^^ nonhelical MHD, 
only n^^is effective, while in helical MHD both kinetic and magnetic helicities play an important role in the growth 
of ME. In the current kinematic models of planetary magnetism [J , magnetic field is generated by kinetic helicity, 
for which planetary rotation (spin) plays an important role. These models appear to work for all the planets except 
Mercury, which rotates far too slowly. We conjecture that 11^^ (independent of helicities) probably plays an important 
role in the generation of magnetic field of Mercury. 

In helical MHD. the helical contribution to the magnetic energy growth goes as [see Eq. (|49|) ] 



dE'' 
IT 



aril - brMTK (58) 



where a and h are positive constants. The term ar\i is always positive independent of the sign of Hm, but —brMfK is 
positive only when Hm and Hk are of the opposite sign. In numerical simulation of Brandenburg |20| and EDQNM 
calculation of Pouquet et al. Q, HmHk < for small k, and H^Hk > for large k. Hence, —brMrx term is positive 
for small k, resulting in positive dE'' /dt. Let us compare the above result with the dynamical dynamo of Pouquet et 
al. 0, Field et al. 0, and Cho 8\. 

The kinematic dynamo predicts that the growth parameter a is proportional to Hk, i.e, a — au oc (u • V x u). 
The kinetic model was generalized by Pouquet et al. Q , Field et al. [3| , and Cho !8| . In absence of a mean magnetic 
field they find that 

a fv au + ab = ... (u • V x u) — ... (b • V x b) , (59) 

where ... denotes certain time scales (which are always positive). It implies that a gets positive contribution from 
both the terms when Hm and Hk are of opposite signs. This result is consistent with Eq. H58|) . The direct numerical 
simulation of Pouquet and Patterson [l9| indicate that Hm enhances the growth rate of ME considerably, but that is 
not the case with Hk alone. This numerical result is somewhat inconsistent with results of Pouquet et al. and others 
13 (Eq. (jS^ ). but it fits better our formula lf3S|l {dE'^/dt = if tm = 0). Hence, our formula (|^ probably is a better 
model for the dynamically consistent dynamo. 

In the following section we will summarize our results. 

IV. CONCLUSIONS 

In this paper we have applied first-order perturbative field theory to calculate the renormalized viscosity, renormal- 
ized resistivity, and various cascade rates for helical MHD. We find that the renormalized viscosity and resistivity are 
unaffected on introduction of both kinetic and magnetic helicities. Our result is consistent with Zhou's calculation 
|27| for helical fluid turbulence. 

We find that the energy cascade rates get significantly altered by helicity. Since magnetic helicity is a conserved 
quantity in MHD, Frisch |25| had argued for k~^ energy spectra at small wavenumbers. However, in this paper we 
calculate energy fluxes in the Kolmogorov's inertial range, where we find direct energy cascades. The fluxes are shown 
in Table ITU The main results of our calculation are as follows: 

1. The magnetic energy flux has two components: (a) the nonhelical part which is always positive, (b) the helical 
part which is negative (assuming HmHk < 0). The inverse cascade resulting due to helicities is consistent with 
the results of Pouquet et al. |5|, Brandenburg ^^, and others. 

2. The u — u flux HJ^^gets a inverse component due to kinetic helicity. This implies that KE flux decreases in 
presence of Hk , a result consistent with that of Kraichnan [Sfl • 

3. The growth of large-scale magnetic fleld in the initial stage of evolution results from ^^^ +^h'^figiicai~^^b<heHcai ■ 
In this paper we have computed the relative magnitudes of all three contributions, and find that all of them to 
be comparable, although H^^ is somewhat higher. Pouquet et al. [5||, Pouquet and Patterson il9|, Brandenburg 
|20l |. and many others highlight H^J^^^;,-^^; transfer, and generally do not consider H^_5and 11^^^^^^^^^ fluxes. 

4. Regarding positive Hm, the flux of magnetic helicity TIhm is backward. 

Most of the earlier papers (e.g. Pouquet et al. |a]) assume Alfven time-scale to be the dominant time-scale for MHD 
turbulence. We have taken nonlinear time-scale (based on Kolmogorov's spectrum) to be the relevant time-scale based 
on recent numerical [IE [H [ill and theoretical [111 13 E EUlll work. 
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Using the flux results we have constructed a simple nonlinear and dynamically consistent galactic dynamo. Our 
model shows an exponential growth of magnetic energy in the early phase (much before saturation). The growth 
time-scale is of the order of 3 x 10® years, which is consistent with the current estimate j^. In our paper we have 
discussed the growth of magnetic field at scale comparable to forcing scales. In real dynamo, the magnetic field at 
even larger scales also grow. This growth may be due to inverse cascade of magnetic energy. This problem is beyond 
the scope of our paper. 

Some of the results presented here are general, and they are expected to hold in solar and planetary dynamo. For 
example, we find that LS velocity field supply energy to LS seed magnetic field. This is one of the sources of dynamo. 
Hence, if we solve first few MHD modes in spherical coordinate with kinetic forcing, it may be possible to capture 
some of the salient features of solar and planetary dynamo. 

In summary, the energy flux studies of helical MHD provide us with many important insights into the problem of 
magnetic energy growth. Its application to galactic dynamo yields very interesting results. A generalization of the 
formalism presented here to spherical geometry may provide us with insights into the magnetic field generation in the 
Sun and Earth. 
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APPENDIX A: CALCULATION OF (S^^(k' p q)) 

The expressions for (^S^^') for helical MHD are 

+n{k,p, q)G''^{p, t - t')C"(fc, i, t')C"(<z, i, t') 

-,T^(.,.,.)G-(.,t-0^^4^^^^ 

+T9(fc,p, g)G""(g, t - t')C^{k, t, t')C^{p, t, t') 

^rj.,/, .r^uui , ,,.HK{k,t,t')HKip,t,t') 
+Tg{k,p,q)G {q,t-t) 



^2 p2 

{S-\k'\p\q)) = I dt'{2nf[T2{k,p,q)G-"{k,t-t')C\p,t,t')C\q,t,t') 
+T^{k,p,q)G^^{k,t ~ t')HM{p,t,t')HM{q,t,t') 

+Tr{k,p, q)G''\p, t - t')G^{k, t, t')G\q, t, t') 

+n{k,p,q)G'\p,t-t' f''^^^y^ HM{q.t,t') 
+T^^{k,p,q)G^^{q,t-t')G^{k,t,t')C\p,t,t') 

+Ti,ik,p,q)G-^iq,t-t')^^^^^^^HM{p,t,t') 
{S'-ik'\p\q)) = f dt'{27rf[Ts{k,p,q)G'\k,t~t')G-{p,t,t')G'{q,t,t') 



(Al) 



(A2) 



..HK{p,t,t') 



+n{k,p,q)G''\k,t-t') "^^; ' ' HM{q,t,t') 

pZ 

+Te{k,p, q)G''^{p, t - t'jG'ik, t, t')G\q, t, t') 

+T^ik,p,q)G^''{p,t-t')HM{k,t,t')HM{q,t,t') 
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+Ti2{k,p,q)G'"'{q,t ~ t')C\k,t,t')C^{p,t,t') 



{S'\k'\p\q)) = I dt'{27rf[Uk,p,q)G'\k,t~t')C\p,t,t')C^{q,t,t') 

,.HK{q,t,t') 



+Ti{k,p,q)G"'{k,t ^ t')HM{p,t,t') 

+Tsik,p, q)G''\p, t - t')G\k, t, t')C''{q, t, t') 

+T^{k,p,q)G'\p,t^t')HM{k,t,t')^^^^^^ 

+Tio(k,p, q)G^''{q, t - t')C\k, t, t')C\p, t, t') 
+Tioik,p,q)G''^iq,t-t')HM{k,t,t')HMip,t,t')] 



(A3) 



(A4) 



where Ti{k,p, q) are given in paper I. To obtain T,-{k,p, q) (helical part) we replace all the second rank tensors of the 
type Pjaik) by ejaih. 

APPENDIX B: CALCULATION OF (S^^^(k'|p|q)) 

The quantity (S'^'^(k'lplq)) of Eq. ^ simplifies to 



(5^-(k'|p|q)) = -^h,m{b^{k')u,ip)b^iq)) 



T 



-^jlm-hr {Ui{q)hrn{k')bj{p)) 



+<^jlm.-^ {bi{q)brn{k')Uj{p)) 

which is computed perturbatively to the first order. The corresponding Feynman diagrams are 



{S"^'{k'\p\q)) 




^C^-r>^./'>^y\y\j^ 




+ 



(Bl) 








(B2) 



Here empty, shaded, and filled triangles (vertices) represent eijm,—£ijmkiki/k'^ and eij^hki/k'^ respectively. The 
empty and filled circles (vertices) denote (— i/2)Pj~,^ and —iP^,n respectively. The solid, dashed, wiggly (photon). 
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and curly (gluons) lines denote {uiUj),{bibj),G'^^, and G^^ respectively. When we substitute {uiUj) , {bibj) using 
Eqs. (|18I15|I . we obtain terms involving C^{p,t,t')H(^M.K){<l:'tTt')- The resulting expression for (5^*'(k'|p|q)) is 






(5^-(k'|p|q)) = / dt' (27:) 

J — OO 

+T32 (fc, p, q)G'''' (fc, t - t')C^^ {p, t - t')HM (g, t - t') 
+T33ik,p, q)G"''{p, t - t')HM{k, t - t')C''\q, t - t') 
+T34 (fc, p, g)G™ (p, t-t')C''\k,t- t')HM {q, t - t') 
+n5{k,p,q)G''\q,t- t')HM{k,t-t')C^'^{p,t-t') 

..HK{p,t-t') 



n2 



+T,e{k,p,q)G'\q,t-t')C''\k,t-t') 
+T37{k,p,q)G''Hk,t~ t')HM{p,t~ t')G'''^{q,t~ t') 
+ns{k,p, qjG'Hk, t - t')C'"'{p, t - t') ^^^'^'f ~ ^'^ 
+T39ik,p,q)G'\p,t^ t')HM{k,t- t')C^-^{q,t- t') 
+T^oik,p, q)G''\p, t - t')C'"{k, t - t') HK{q,t-t') 

q 

+Tiiik,p, 9)G""(g, t - t')HMik, t - t')G''\p, t - t') 
+T42(fc,p, 9)G™(9, t - t')G''\k,t- t')HMip, t ~ t')] 

+Tisik,p, q)G''\k, t - t' f''^P'l'^'^ C"'{q, t - t') 

pZ 

+n^{k,p, q)G''Hk, t - t')C^"{p, t - t')HM{q, t - t')} 
+r45(fc,p, 9)G™(p, t - t')HM{k, t - t')C'^\q, t - t') 

+Ti^(k, p, q)G^^ (p, t-t')C'\k,t- t')HM (g, t - t') 
+T^7{k,p,q)G''\q,t ^ t')HMik,t - t')G^^{p,t ^ t') 

+T48(fc,p, q)G'\q, t - t')C''\k, t - t' fK{p,t-t') 

pz 



(B3) 



The terms T'3i..48(fc,p, q) can be obtained in terms of antisymmetric tensors ejai, P^m ^^'^- They have not been listed 
here due to lack of space. 
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TABLE I: The values of I^ = (n^/n)/(A'")^-^ calculated using Eqs. II47I48II for Alfven ratios ta = 1 and ta = 5000. 



rA 



TA = 5000 



TU< 
^U> 
JU< 

^b> 

Tb< 
^U> 
Tb< 

^b> 

TU< 

^b< 



0.18- 

0.54 



0.19 - O.lOrf, 



2.04rl, + l.QSrKTM 
- 1.9rli + 2.02rKrM 



-25rM + 0.35rK 



0.53 - 0.28r| 



1.9 X lO""* + 1.4 X 10~Vif + 2.1 X IQ-Vxr-M 
-5.6 X 10-"^ - 1.1 X 10-Vif + 5.4 X lO-^r^rM 
1.4 X 10~* - 1.02 X lO^VIf + 5.4 X lO^^'r^rM 

-4.1 X lO-^Af + 8.1 X IQ-^rK 



TABLE IL The values of energy flux ratios Ily/n for various values of rA,rK, and tm for k ^" region. The first and second 
entries are nonhelical and helical contributions respectively. 



{rA,rK,rM) 



(5000,0.1,-0.1) 

(5000,0.1,0.1) 

(1,0.1,-0.1) 

(1,0.1,0.1) 

(1,1,-1) 

(1,1,1) 

(1,0,1) 



n::$/n 



^b> 



/n 



Kt/n 



(2.5 X 10-"*, 

-9.4 X lO"'*) 

(2.5 X 10"*, 

9.4 X 10-*^) 

(0.37, -0.027) 

(0.35,8.3 X 10"'' 

(0.083, -0.60) 

(0.34,0.081) 

(0.23, -0.8) 



(1.0,-0.0053) 

(1.0,-0.0053) 

(0.13,-6.9 X 10"' 

(0.12,-6.5 X 10"' 

(0.029, -0.015) 

(0.12, -0.064) 

(0.081,0) 



(3.2 X 10"*, 

-3.7 X lO"'^) 

(3.2 X 10"*, 

3.7 X lO"'') 

) (0.43,-4.4 X 10"*) 

) (0.4,8.1 X 10"*) 

(0.095,-9.9 X 10"^) 

(0.39,0.079) 

(0.26,0.013) 



(-9.7 X 10-^ 
-9.0 X 10"^) 
(-9.7 X 10"^ 

9.0 X 10"^) 
(0.13,-0.027) 
(0.12,-7.7 X 10" 
(0.028, -0.61) 
(0.12,-0.075) 
(0.078, -0.86) 



